Astronomy & Astrophysics manuscript no. rxj09 


©ESO 2013 


March 28, 2013 





Resolving the molecular gas around the lensed quasar 

RXJ091 1.4+0551* 

P. T. Arm 1 ' 2 ' 3 ' 4 , F. Boone 1 ' 2 , D. T. Hoar 3 ' 4 , P. T. Nhung 3 , A. WeiB 5 , J.-P. Kneib 6 , A. Beelen 7 , and P. Salome 8 

1 Universite de Toulouse; UPS-OMP; IRAP; Toulouse, France 

2 CNRS; IRAP; 9 Av. colonel Roche, BP 44346, F-31028 Toulouse cedex 4, France 

3 VATLY, INST, 179 Hoang Quoc Viet, Cau Giay, Hanoi, Vietnam 

4 Institute of Physics, 10 Dao Tan, Hanoi, Vietnam 

5 Max-Planck Institut fur Radioastronomie, Auf dem Hugel 69, 53121 Bonn, Germany 

6 Laboratoire d'Astrophysique de Marseille, CNRS -Universite Aix-Marseille 

7 Institut d'Astrophysique Spatiale, Bat. 121, Universite Paris-Sud, 91405 Orsay Cedex, France 

8 LERMA, Observatoire de Paris, UMR 81 12 du CNRS, 75014, Paris, France 

Received ; accepted 

ABSTRACT 

We report on high angular resolution observations of the CO(7-6) line and millimeter continuum in the host galaxy of the gravitation- 
ally lensed (z~2.8) quasar RXJ09 11.4+0551 using the Plateau de Bure Interferometer. Our CO observations resolve the molecular 
disk of the source. Using a lens model based on HST observations, we fit source models to the observed visibilities. We estimate 
a molecular disk radius of l+0.2kpc and an inclination of 69+6 deg, the continuum is more compact and is only marginally re- 
solved by our observations. The relatively low molecular gas mass, M gas = (2.3 + 0.5) x 10 9 M©, and far-infrared (FIR) luminosity, 
Lfir = (7.2 + 1.5) x 10 11 L©, of this quasar could be explained by its relatively low dynamical mass, M dyn = (3.9 + 0.9) x 10 9 M©. 
It would be a scaled-down version the quasi-stellar objects (QSOs) usually found at high-z. The FIR and CO luminosities lie on the 
correlation found for QSOs from low to high redshifts and the gas-to-dust ratio (45 + 17) is similar to the one measured in the z=6.4 
QSO, SDSS Jl 148+5251. Differential magnification affects the continuum-to-line luminosity ratio, the line profile, and possibly the 
spectral energy distribution. 



1. Introduction 

In the past 15 years, submillimeter/millimeter (submm) observa- 
tions of the interstellar medium in high-redshift (high-z) galax- 
ies have opened a new window for studying the early Universe 
and the evolution of galaxies. In particular, continuum emission 
from dust heated by intense star formation (> lOOM yr _1 ) as 
well as large amounts of molecular gas (> 10 10 M©) have been 
observed in quasi-stell ar objects (QSOs) back to the end of the 
epoch of reionization ( [Solomon & Vanden Bout|[2005| [Carilli 



|& W alter 2013). Simultaneously, a populati on of high-z galax- 
ies dubbed submillimeter galaxies (SM Gs, Blain et aT| 2002 ; 
|Ivison et al.|201 1| Riecher s'et al.|201 fb] ) with similar molecular 
gas masses, star formation rates (SFR) and evolution have been 
discovered. Recent observations suggest a possible evolutionary 
link between these two populations. The large gas reservoirs ob- 
served in SMGs are supposed to be driven to the galaxy center, 
feeding intense starbursts and leading to the growth and activity 
of the central black holes that would become QSOs. The SMGs 
would thus transform into QSOs on a time-scale of -lOOMy r 
( [Carilli & Wang|2006[|Coppin et al.|20"08l|Simpson et al.|2012| ). 
If confirmed, this scenario would link the history of cosmic star 
formation with that of black hole growth, and it would give new 
insights into the role of feedback mechanisms in galaxy evolu- 
tion. However, more detailed studies of individual galaxies are 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



required to verify whether the morphologies and star-forming 
properties are consistent with this scenario. 

Gravitational lensing is an essential tool for conducting such 
studies at sensitivities and physical resolutions that would other- 
wise be out of reach. The spatial stretching induced by lensing 
combined with high-resolution instruments allows us to obtain 
information on the morphology and the kinematics of the source. 
The pi oneer obse rvation of the Cloverleaf using the CO(7-6) 
line ( All oin et al.||1997| ) was the first to demonstrate the power 
of such studies. Resolving the source is also important for cor- 
rectly interpreting the integrated fluxes of lensed systems be- 
cause the regions observed in line and continuum or at different 
frequencies may not overlap and may therefore undergo differ- 
ential magnification. 

In this paper we present high-resolution Plateau de Bure 
Interferometer (PdBI) observations of the CO(7-6) line and the 
millimeter continuum emission in the quadruply lensed quasar 
RXJ091 1.4+0551 at z = 2.796 ([Bade et al.|1995| |Hagen et al.| 
T995j|Bade et"aL][T997l |WeiB et al.||2012] TThis QSO is partic- 
ularly interesting because the lens, which includes the contribu- 
tion of an intervening galaxy cluster, has previously been well- 
studied ( [Burud et aL]|1998[ |Kneib et all|2000| ), and recent ob- 
servations revealed a molecular gas mass M gas = 4.2 x 10 9 M© 
(Hainl ine et al. 2004[) and a far -infrared (FIR) luminosity L FIR = 
0.8 x 10 1Z L G |wu et alT]|2009| ) that are an order of magnitude 
lower than usually observed in QSOs. By measuring the source 
size and morphology we aim at measuring the dynamical mass 
to find out whether these low values are measured because this 
is a lower mass galaxy or because this is an evolved galaxy that 
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Fig. 2. Left panel: observed velocity gradient. The grayscale 
image is a zoom into the CO(7-6) map shown in Fig.[T] The 
blue and red contours correspond to the CO(7-6) line emis- 
sion integrated in the ranges [-135, -15] and [-15, 105] km s" 1 , 
respectively. The contours start at 3-cr and are spaced by 
cr=104mJybeam -1 kms -1 . 2nd to 4th panels: Simulated ve- 
locity gradients with two lensed circular Gaussian sources of 
FWHM p= 120 mas, separated by 110 mas and with a position 
angle o f 10 8, 18 and 63 deg, east from north, respectively (see 
SectionlO). 




Fig. 3. Four quasar images obtained with a lens model built fol- 
lowing |KneiF£Far](|2^0]). The HST image is overlaid in ma- 
genta contours. The red lines represent the critical lines and the 
two yellow circles show the positions of the two lensing galax- 
ies. 



has consumed its gas already. We also investigate the effects of 
differential magnification on the continuum/line luminosity ratio 
and on the spectra. We adopt the ACDM concordance cosmol- 
ogy: H = 71 kms" 1 Mpc" 1 , D M = 0.27 and Q A = 0.73. 

2. Observations and results 

2.1. Observations 

The data were collected with PdBI during ~18 hours on-source 
in good observing conditions in February and March 2007. The 
receivers were tuned at 212.485 GHz. Pointing and focus mea- 
surements were performed at intervals shorter than one hour. 
Several quasars within 10 deg from the source were used as 
calibrators: MCW349, 0851+202, 0906+015, 0923+392, and 
J04 18+380. Data reduction was performed using the GILDAS 
standard package ( jhttp : //www . iram . f r /IRAMFR/ GILDAS/] ). 
The data were averaged in 46 velocity bins of 30 kms" 1 each. 
The dirty map was reconstructed on a square grid with a pixel 
size of 80 mas using natural weighting. The map was cleaned 
with the CLEAN method (H6gbo m|1974| ). The clean beam is an 
elliptic Gaussian with a full width at half maximum (FWHM) of 
0.72" and 0.27" along the major and minor axes, respectively, 
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Fig. 5. Upper panel: CO(7-6) map shown in the Fig.[l] Middle 
panel: The best- fit Gaussian circular source with p= 151 mas 
lensed and convolved by the PdBI clean beam. Lower panel: 
The best-fit elliptical source lensed and convolved with the 
PdBI clean beam. The ellipse parameters are p ma j=260mas, 
Pmin/Pmaj=0-36, and #=108 deg. The left insets show the resid- 
uals at the A images, the contours are at lx<x and 2x<x. 




Fig. 6. Color image representing the best-fit elliptical Gaussian 
source to the CO(7-6) emission. The yellow curve shows the 
caustic line derived from the lens model. To quantify the effect 
of differential magnification we considered three virtual sources. 
SI, S2, and S3 are magnified by 18.8, 39.3, and 13.8, respec- 
tively. 



and a position angle of 15 deg from north to east. The noise level 
is 1.7 mJy beam -1 in each velocity channel. 

2.2. Flux, morphology, and velocity gradient 

The integrated CO(7-6) line map, the continuum map, and the 
integrated spectrum are shown in the Figure[T] The HST observa- 
tions of the quasar with the astrometry calibrated to match that 
given by Kneib et aL] ( |2000| ) are overlaid in contours. The line 
width, FWHM=120+14kms-\ the line peak, S max =48+ 5 mJy 
and the line redshift are consistent with those obtained by WeiB 
et al. (2012 ) from independent PdBI observations. However, due 
to a higher noise level we do not detect the weak line wings 
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Fig. 1. Upper left: CO(7-6) emission integrated over the velocity range [-120, 60] kms -1 . The black contour levels start at 3x<x 
and are spaced by cr = 135 mJy beam -1 kms -1 . The white ellipse at the bottom right represents the beam at half maximum 
(0.71"x0.28"). The HST map contours are overlaid in cyan. Upper right: The 212 GHz continuum map. The black contour levels 
start at 3x<x and are spaced by cr = 0.28 mJy beam -1 . Lower left: The CO(7-6) integrated spectrum with the best Gaussian fit 
overlaid. Lower right: best-fit mo dified black body spectral energy distribution to the FIR continuum measurements derived by Wu 
|eTaLl ( [2009l ), |Barvainis & Ivison| ( |2002] ), and us (red dot at A =1.41 mm). 
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Fig. 4.x 2 variations with the source parameters. The upper panels show the variations of ;^ 2 -24755 for the line maps and the lower 
panels show ^ 2 -22 188 for the continuum maps. From left to right the abscissae are the source size (FWHM), p, assuming a circular 
Gaussian source; the offset in RA, x ff; the offset in Dec, y ff ; the position angle 6 assuming an elliptical Gaussian source; the source 
major axis, p ma j, and the source axis ratio r = p m in/Pmaj- The horizontal line shows the 68% confidence level (1 x cr) taking into 
account the number of free parameters. 



reported by these authors that could be the signature of an out- 
flow. The line intensity is 7co(7-6) = 6.1 ± 1.0 Jy kms -1 , which 
translates into a line luminosity ^0(7-6) = (4-36 ± 0.8) x p~ l x 
10 10 Kkms" 1 pc 2 , where \i is the lensing magnification factor. 
The continuum flux, S212GHZ ~ (7-4 ± 0-9) x p~ l mJy, is also 
consistent with other submm continuum measurements as shown 
Fig.[T] The best-fit modified black body spectral energy distribu- 



tion (SED) with ft = 1.5 gives a dust temperature r d = 39 ± 3 K 
and a total FIR luminosity L FIR = 273 ± 45 x pT l x 10 11 L . 

The line and continuum emissions show similar morpholo- 
gies that are consistent with the configuration of the four quasar 
images observed with HST. The B image is tentatively detected 
at the 2.6-<T level, consistently with the B/A flux ratio of 0.21 
measured by Wei B et al. ( 2012| ). The overlay also suggests that 



3 



P. T. Anh et al.: Resolving the molecular gas around the lensed quasar RXJ091 1.4+0551 



the PdBI observations are offset with re spect to HST observa- 
tions, as was noted by WeiB et al. (2012). 

Fitting four point sources (i.e., beams) with relative positions 
fixed to those observed with HST to the CO emission leads to 
large residuals. This suggests that the molecular gas around the 
quasar is resolved. This hint is further supported by the evidence 
of a velocity gradient across the map. Indeed, as shown in the 
Fig|2j the emission corresponding to the blue half of the line is 
spatially offset from the emission corresponding to the red half 
of the line. The offset between the peaks of the blue and red 
emission is ~0.3" in the northwest direction, which is close to 
the beam width in this direction (0.28"). The offset is therefore 
significant and suggests a resolved dynamical structure. 

2.3. Source shape modeling 

To interpret the PdBI observations in terms of source proper- 
ties we have built a lens model using the LENSTQQL software 
flKneib et al.|1996l| Julio et al.|2007H Julio & Kneib|2009) follow- 
ing the same approach as |Kneib eTa l. ( 2000). The model consists 
of three elliptical isothermal sphere potentials corresponding to 
the nearby galaxy cluster, the main lensing galaxy, and a smaller 
neighboring galaxy. The four quasar images obtained are shown 
in the Fig. [3] with the HST observations overlaid in contours. 
In parallel and as a consistency check we compared our results 
with those obtained using a much simpler gravitational potential 
dBurud et al.|1998| ) and a different software ( |Hoai et al.|2013| ). 

To derive the source properties we proceeded by model- 
fitting in the Fourier plane. This method allows for a rigorous 
treatment of the errors and avoids artifacts inherent to imag- 
ing with interferometer data (see e.g. Bussm ann et al. 2012 ^. 
For each source parameter value the sky image produced by 
LENSTOOL is Fourier-transformed and sampled at the (u, v) co- 
ordinates of the visibilites measured by PdBI. The source param- 
eters are varied to minimize^ 2 = £f ||V° bs - V^ od \\ 2 / erj , where 
N = 13710 is the number of visibilities measured by PdBI, V° bs 
and V™ od are the measured and modeled complex visibilities, 
and er i are the noise root mean squares (RMS) of the measure- 
ments. In addition to the source shape parameters we considered 
an astrometric offset in RA and Dec (x ff and y fi) and a scaling 
factor (/, equivalent to a source flux parameter). Any significant 
offset of the CO and/or the continuum emission from the QSO in 
the source plane was excluded because it would produce a dif- 
ferent configuration for the images. Such an offset is generally 
observed in merging systems, which does not seem to be the case 
here given the single narrow line. The number of free parame- 
ters was taken into account when computing the uncertainties 
following [PFess et al.| | |2002) . 

We first used a simple circular Gaussian source model for 
which the only shape parameter is p, the FWHM of the Gaussian. 
The best-fit values obtained arep = 151+20 and 40 + 8 mas for 
the line and continuum emission, respectively. The offsets ob- 
tained are x off = 0.15 + 0.04 and y oS = 0.22 + 0.08", consistent 
with the uncertainties in the astrometry of the HST map. The 
corresponding x 2 curves and best-fit sky images are shown in 
the Figs.|4]and[5] The offsets obtained in the other fits described 
below are similar within the uncertainties. These results show 
that the region emitting the CO(7-6) line is clearly resolved and 
that the region emitting the continuum is resolved with a confi- 
dence of ~ 5er. 

A galaxy is not expected to be circular on the sky. However, 
fitting a non-circular model would require fitting two more pa- 
rameters (elongation and orientation), which seems too much for 



the signal-to-noise ratio and the resolution of our data. We there- 
fore followed a two-step approach. We first fitted the orientation 
of an elliptical Gaussian assuming a fixed and arbitrary elonga- 
tion for the source, namely an elliptical Gaussian source with 
an axis ratio of p m in/Pmaj = 0-5. To keep the size of the source 
consistent with p, the FWHM obtained with the circular source 
fit, we derived pmin and p maj from (p 2 ^ + p^ aj /2) 1/2 = p. We 
thus obtained a position angle of 6 = 108 + 20 and 104 + 71 deg 
from north to east for the line and the continuum, respectively 
(Fig. [4]). It is remarkable that for the line emission the minimum 
X 2 reached is significantly lower than for the circular source best- 
fit ( A;^ 2 =10, i.e., a sig nificance of 95.5% for four free parame- 
ters, [Presretar[2002]), indicating that the source is indeed elon- 
gated. This is not the case for the continuum, and the large uncer- 
tainty on 6 shows that the position angle of the continuum source 
cannot be constrained by our data. We checked that changing the 
source axis ratio does not affect the value of 6 obtained for the 
line emission within the uncertainties. 

In a second step we fixed the position angle to 6 = 108 deg 
and we simultaneously fitted the major axis FWHM of the el- 
liptical Gaussian, p ma j and its axis ratio, r = Pmin/Pmaj- This 
was done for the line-integrated map only because the contin- 
uum does not appear to be elongated. The results are p ma j = 
260 + 40 mas and r = 0.36 + 0.09, excluding a circular source 
with a confidence > ler. The improvement with respect to the 
circular source is also noticeable in the sky image (Fig. [5]). We 
checked that these results are robust against small variations of 
the lens model parameters within the uncertainties. 

If this is a rotating disk, there should be a velocity gradi- 
ent along the major axis. To test whether this solu tion is con- 
sistent with the velocity gradient noticed in Section |2^2{ we re- 
placed the elliptical model by two circular Gaussian sources of 
FWHM, p = 120 mas, separated by 110 mas along the major 
axis and assumed that the northern source is red and the other is 
blue. The result is shown in the second panel of Fig. [2] As can 
be seen, the velocity gradient thus produced looks similar to the 
observed one. This would not be the case for other orientations, 
as illustrated in the two other panels. This good agreement with 
the kinematics information supports our fit results. 



3. Discussion 

From our best-fit models we derived global magnification fac 
tors of \i- 18. 1+3.6 and 37.7+7.2 for the line and the continuum 
respectively, yielding ^0(7-6) ~ (2.42 
andS 2 i2GHz = (196 + 30)/Jy. 



0.50)xl0 9 Kkms" 1 pc 2 



Assuming yu=22, |Riechers et aTT| ( [201 laj ) measured ^co(i-o) = 
3.39+0.48 Kkm s _1 pc z , which implies ^co(i-o)/ L co(7-6) = L2if 
we assume yu=18.1 for both CO lines. These authors also noted 
that the corresponding gas mass M gas = (2.3 + 0.5) x 10 9 M Q 



(usin g the conversion factor a = 0.8 from |Downes & Solomon| 
1998 ) is at the low end of the observed range of gas masses in 
QSOs. 

Assuming that the submm continuum measurements are all 
emitted from the same region (no differential magnification) and 
that they can be modeled with a single modified black body SED 
as in Fig.[T] the total FIR luminosity is L FIR = (7.2 + 1.5) x 
10 11 L n . If the FIR is entirely prod uced by star forma tion, it cor- 



responds to an SFR= 1 10 M Q yr 1 ( [Kennicutt 



1998b a). From the 



best- fit SED and making the same assumptions as |Boone et al.| 
( [20071 ), we derive a dust mass M d =(5.9 + 1.2) x 10' M Q , corre- 
sponding to a gas-to-dust ratio M gas /Md = 45 + 17, i.e., similar 
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to that me asured in the z=6A QSO SDSS Jl 148+5251 (Valiante 
|et al.|2012» . 

The FIR and CO luminosities are compatible with the L FIR - 
Lp n correlation me asured in quasars from low to high redshifts 
by Riechers (2011) and do not match the correlation found by 
Bothwell et al.| ( |2013| ) for SMGs, L' co being a factor ~ 3 too 
low. KXJ (FTTs host properties would therefore be consistent 
with the transition scenario from SMGs to QSOs. 

Ignoring the differential magnification, i.e., using the same 
magnification factor for both the CO and the FIR emission, 
would have lead to an overestimate of the star formation effi- 
ciency measured by L FIR /L^ . This case clearly illustrates why 
this effect needs to be taken into account, in particular in QSOs 
where the FIR may come from a more compact region than the 
CO lines. The star formation efficiency may thus have been over- 
estimated in several unresolved lensed QSO hosts. 

The best-fit source shape parameters can be used to ob- 
tain a tentative estimate of the geometry. At z = 2.796, 1" cor- 
responds to 7.98 kpc, the characteristic radius of the CO disk 
is therefore R = 1 + 0.2 kpc, and its inclination angle is i = 
90 - arcsin(r) = 69 + 6deg. This high inclination seems to 
contradict the remarkably narrow line observed, which suggests 
a face-on disk. However, the corresponding dynamical mass, 



M dyn = R x Av 2 x G 1 x sin i 2 



(3.9 + 0.9) x 1O 9 M , is a 



factor 2 larger than the molecular gas mass, implying a gas frac- 
tion of 50%, similar to the one observed in other high-z QSOs 
(e.g., |Riechers et al.|[2009l |Solomon & Vanden Boutff2005] ). In 



other words, the line may be truly narrow (the line width cor 
rected for inclination would be W = 129 + 20kms _1 ) because 
this quasar host may lie on the 'light' end of the mass distribu- 
tion of quasar hosts ( [Caril li & Wang 2006). This interpretation is 
also consistent with the X-ray luminosity, Lx = 1 .4x 10 44 erg s" 1 
pan et al.||2009 assuming the same magnification factor as 
for the FIR, i.e., 37.7) because according to the well-known 
^bh - cr+ relation in local elliptical galaxies (Gultekin et al. 



2009| ) and following |Bothwell et al.| ( |2013| ), the line width ex- 
pected is W ~ ov/1.5 " 135kms- A for M BH = 1O 7 - 7 M . Thus, 
RXJ091 1 appears to be a scaled-down version of the QSOs usu- 
ally found at high-z, residing in a lighter halo, containing less 
gas, forming fewer stars and having a lighter central black hole. 

However, the line may still be broader than observed due to 
the effec t of differen tial magnification. Indeed, as was recently 
noted by |Hoai et al.|p013] ), the source is near the caustic cusp 
and differential magnification is expected to be significant. In 
fact, as shown in Fig. 6] the caustic lines cross our best-fit source 
model mainly along the minor axis where the line-of- sight veloc- 
ities should be close to the systemic velocity. Differential mag- 
nification should therefore boost the center of the spectral line 
when integrated over the entire source, resulting in a narrower 
profile. To illustrate this effect we computed the magnification 
factor of three virtual circular Gaussian sources separated by 
100 mas and located along the major axis, as shown in Fig. [5] 
We find that the central source S2 is more magnified than SI 
and S3, by a factor 2.1 and 2.8, respectively. Correcting the line 
profile for this effect would require estimating the rotation curve 
from higher resolution observations. A higher sensitivity is also 
re quired to investigat e whether the broad line wings discovered 
by |WeiB et al.| (2012 ) could be caused by this effect. In a similar 
study of another lensed quasar, |Riechers et al.| ( |2008| ) found a 
much weaker differential effect on the CO line. This is because 
the CO region observed in their source is much more extended 
than the region delimited by the caustic lines; this is not the case 
inRXJ0911. 



This case study shows how gravitational lensing can provide 
insights into sources that are scaled-down versions of the usually 
selected ones. But it also underlines the possible bias introduced 
by lensing on the integrated measurements due to differential 
magnification. High-resolution observations are required to cor- 
rect for this effect and to access the physical properties of the 
sources. ALMA will probably allow us to study the molecular 
gas morphology and kinematics in this and other lensed QSOs 
and SMGs in more detail. In particular, ALMA will be essential 
for studying the scaling relations toward the lower intrinsic gas 
masses. Indeed, the few other lensed quasars with similar low 
molecular gas masses lack good lens models (Riechers 2011, 
and references therein). This will allow us to better understand 
the relation between molecular gas mass, star formation rate, and 
total mass, and to identify possible differences between the prop- 
erties of QSO hosts and SMGs in a more robust way. 
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